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The Lande g factor describes the response of an atomic energy level to an external perturba- 
tion by a uniform and constant magnetic field. In the case of many-electron systems, the leading 
term is given by the interaction /j,b (L + 2 S) ■ B, where L and S are the orbital and spin an- 
gular momentum operators, respectively, summed over all electrons. For helium, a long-standing 
experimental-theoretical discrepancy for P states motivates a reevaluation of the higher-order terms 
which follow from relativistic quantum theory and quantum electrodynamics (QED). The tensor 
structure of relativistic corrections involves scalar, vector, and symmetric and anti-symmetric ten- 
sor components. We perform a tensorial reduction of these operators in a Cartesian basis, using 
an approach which allows us to separate the internal atomic from the external degrees of freedom 
(magnetic field) right from the start of the calculation. The evaluation proceeds in a Cartesian 
basis of helium eigenstates, using a weighted sum over the magnetic projections. For the relativis- 
tic corrections, this leads to a verification of previous results obtained using the Wigner-Eckhart 
theorem. The same method, applied to the radiative correction (Bethe logarithm term) leads to a 
spin-dependent correction which is different for singlet versus triplet P states. Theoretical predic- 
tions are given for singlet and triplet 2P and triplet 3P states and compared to experimental results 
where available. 
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I. INTRODUCTION AND OVERVIEW 
A. Few-electron systems and g factor 

The quantum electrodynamic (QED) theory of 
bound systems describes, among other things, 
three "fundamental" characteristic effects that in- 
volve the spectrum of bound systems, namely, 
(i) the Lamb shift, which is the energy shift of 
bound states due to the self-interaction of the elec- 
trons, and due to tiny corrections to the Coulomb 
force law at small distances, (ii) the g factor of 
bound states, which describes the energy shift 
of a bound state due to the interaction with an 
external, uniform magnetic field (Zeeman effect), 
and (iii) the hyperfine splitting, which is given by 
the interaction of bound electrons with the nu- 
clear magnetic moment. These effects seem to be 
the three most commonly studied QED effects for 
bound states, because of prominent high-precision 
experiments in all three mentioned areas. The 
leading QED corrections to all three mentioned ef- 
fects are given by the self-energy of the orbiting- 
particle, and by vacuum polarization. 

The long standing discrepancy between theory 
and experiment for the Zeeman coupling factor g' L 
for the 2 3 P state of helium [1] has motivated a 
number of independent theoretical papers [2-5] on 



the subject. Here, by convention, g' L is the com- 
plete orbital part of the gj factor for the helium 
P state, including relativistic and radiative correc- 
tions. For hydrogenlike systems, the self-energy 
corrections to the gj factor and to the hyperfine 
splitting can be formulated in a similar frame- 
work [6, 7], by observing that they can be described 
as a "dressed" self-energy correction in an addi- 
tional magnetic field, namely, for the case of the gj 
factor, in a uniform external magnetic field, and, 
for the case of the hyperfine splitting, in the mag- 
netic dipole field of the atomic nucleus. 

For more complex atoms and ions, the theory 
of the g factor is more complicated because in 
higher order, the electron-electron interaction is 
intertwined with the coupling to the external mag- 
netic field. In leading order, the total orbital angu- 
lar momentum L and the spin angular momentum 
S couple to the external magnetic field B as de- 
scribed by the Hamiltonian matrix element 



(H 



^ B (L + 2S)-B^=gj^ B Bfi, (1) 
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where gj is the Lande g factor, and ijlb is the Bohr 
magneton, that is, /is = — e/(2m) where m is the 
electron mass and e = — |e| is the electron charge. 
The orbital angular momentum L and the spin an- 
gular momentum S are summed over all electrons. 
As long as the separation into terms proportional 
to L ■ B and S ■ B remains valid, this gives rise to 
an orbital jj, w 1 factor and a spin g$ ~ 2 factor, 



2 



so that 



.9.7 = 9L 



J(J+1)+L(L + 1)- 5(5 + 1) 



5s 



2J(J+1) 

J(J + l) + 5(5+l)-L(£+l) 
2J(J + 1) 



(2) 



In leading order, the Lande g factor is thus given 
by the well-known formula 

3 J( J + 1) - L(L + 1) + S(S + 1) 
' 9J " 2J(J+1) ' (3) 

In higher order, due to spin-orbit coupling, one 
cannot separate the magnetic-field interaction any 
more into terms proportional to L B and S-B, and 
therefore, one cannot uniquely identify the orbital 
gL and spin gs factors any more. For hydrogen, 
the corresponding mechanism has been discussed 
in Appendix A of Rcf. [8]. The separation into 
gL and gs remains valid up to relative order a 3 , 
where a is the fine-structure constant, provided 
one adds a tiny correction due to a higher-order 
tensor structure, called g x in Refs. [4, 9]. 



B. Angular- momentum algebra 

For P states as opposed to 5 states, the angu- 
lar momentum algebra involved in the calculation 
of the bound-electron g factor can become rather 
complicated, and two approaches have been used. 
In approach (i), used in Refs. [4, 9], the authors 
formulate the entire theory in terms of Wigner 3 J, 
6J, and 9J symbols, which enables them to per- 
form all calculations in terms of reduced matrix 
elements. In turn, these can be written in terms 
of the radial component of the wave functions as 
obtained from variational calculations. 

In approach (ii), which has been used for hydro- 
genlike systems [8] , one first chooses a specific com- 
ponent of the Hamiltonian "vector" hb{L + 25) 
multiplying the magnetic field B, and a specific 
magnetic projection of the reference state. Natural 
choices consist in the z component of the Hamil- 
tonian "vector" and the state with magnetic pro- 
jection [i = | as indicated in Eqs. (15) and (16) of 
Ref. [8] . Due to the Wigner-Eckhart theorem, one 
can then formulate all relative corrections to the g 
factor in terms of ratios, relating the effect calcu- 
lated with a correction to the magnetic Hamilto- 
nian to the leading-order effect, provided one uses 
the same state for each matrix element. This dis- 
entangles the internal atomic degrees of freedom 
from the external degrees of freedom (the mag- 
netic field). For the hyperfine splitting, a similar 
approach is outlined around Eq. (7) of Ref. [7]. 



For helium, it is preferable to formulate the the- 
ory in terms of elements of the radial wave func- 
tions alone, by expressing the matrix elements in 
terms of sums over magnetic projections, where the 
angular and spin degrees of freedom are summed 
over and evaluated in closed form. The latter sum 
can naturally be expressed in terms of a "radial" 
representation of a P state as obtained from a vari- 
ational calculation in a fully correlated, nonrela- 
tivistic basis. Here, we thus choose an approach 
combining ideas from (i) and (ii). First, the rela- 
tivistic and radiative corrections are expressed in 
terms of particular tensor structures, and then, we 
evaluate these on a weighted sum over the pro- 
jections m of the total angular momentum of the 
helium state. This approach combines the advan- 
tages of approach (i), namely, the easy applicabil- 
ity to helium, with the advantages of approach (ii), 
namely, the full disentanglement of the external de- 
grees of freedom (magnetic field) from the internal 
atomic degrees of freedom right from the start of 
the calculation. 

Our investigation is motivated in part by 
an interesting theoretical-experimental disagree- 
ment between the experimental result reported in 
Ref. [1] and theory work described in Refs. [2-4] 
and Sec. V of Ref. [5]. Our calculation is valid 
up to and including relativistic and radiative cor- 
rection of relative order a 3 , and to second order 
in the electron-nucleus mass ratio (for the leading 
nonrelativistic term). We proceed as follows. In 
Sec. II, the terms in the Hamiltonian which govern 
the bound-state g factor are analyzed in terms of 
their tensor structure. The discussion is comple- 
mented in Sec. Ill by an analysis of the spin and 
the tensor reduction of the particular correction 
terms. Finally, in Sec. IV, numerical evaluations 
are described which allow us to obtain a highly ac- 
curate theoretical prediction for the gj factor in 
helium, for 2P and 3P states. Conclusions are re- 
served for Sec. VI. Atomic units with e = h = 1 
(e denotes the physical electron charge) , unit elec- 
tron mass m = 1, a = 1/c and eo = l/(4-7r) are 
used throughout the paper. 



II. HAMILTONIAN 
A. Leading order 

A careful treatment of the g factor requires an 
analysis of the reduced-mass dependence. We de- 
note the electron mass as m and the mass of the 
nucleus as M. The reduced mass /i and the mass 
ratio A are given by 



mM 



M 



X 



' M ' 



(4) 
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The interaction with the external magnetic field, 
in leading order plus the reduced-mass correction, 
is given by [10, 11] 



H M = fJ-B [ql (r a x p„) + y a a ] ■ B 

a 

~ TI ^2(f a x p b ) ■ B . 



(5) 



The finite mass of the nucleus yields a correction 
term (second term) of order 0(A). The sum over a 
and b in Eq. (5) counts the electrons of the bound 
system. The well-known spin factor gs can be 
expressed in the form (including two-loop correc- 
tions) 



9S 



0.328 478 695^- 



(6) 



It is equal to the g factor of the free electron includ- 
ing the anomalous magnetic moment. The terms 
proportional to gi, in Eq. (5) give rise to an or- 
bital factor <7j, = 1 — m/M. The terms in Eq. (5) 
contain all terms of relative order 0(A). We note 
that this scaling of <?l, which has originally been 
derived in Rcf. [10], goes beyond the 'trivial scal- 
ing' of momenta and distances which is discussed 
below. 

Namely, in general, the scaling of the momenta 
and distances with the reduced mass entails the 
scaling factors (see Appendix A) 



(1 + A)- 



(7) 



It results in prcfactors of the form (1 + A)" with 
a certain scaling degree n. For the leading terms 
given in Eq. (5), we have n = 0, and the terms 
commute with the nonrelativistic Hamiltonian of 
the helium atom [without mass polarization; see 
Eq. (A2) below] 
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V 2 ^ 



E- 



(8) 



a>b 



While the first-order correction to the wave func- 
tion due to the magnetic interaction vanishes, the 
mass polarization term i? mp 



mp 



A ^ . . 

r a>b 



(9) 



generates a nonvanishing perturbation to the wave 
function. The perturbation can then be evaluated 
on the leading-order Hamiltonian (5). For the fi- 
nite mass effect of order 0(A 2 ), one additional ef- 
fect is the mass polarization correction to the sec- 
ond term in Eq. (5). The third-order term in- 
volving the leading magnetic interaction term in 
Eq. (5) and two mass polarization insertions also 
yields a finite-mass correction of second order in A. 
Its effect on the g^ prefactor and off-diagonal cor- 
rections are discussed in the following. 



B. Tensor decomposition of the Zeeman 
Hamiltonian 



Let us now turn to the tensor decomposition of 
the Zeeman Hamiltonian. The first term in Eq. (5) 
can be rewritten in the form 



Hmo — fJ-B Go ■ B 7 
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— a er c 



u 0,a 



(f a x p a y , 



do = 1 



(10a) 
(10b) 

(10c) 



where v a is a vector coefficient and do is a diag- 
onal (scalar) coefficient multiplied only by a spin 
matrix. Here and in the following, Cartesian co- 
ordinates arc denoted by superscripts, that is, the 
x component of vo.a is given as Vq a (the super- 
script assumes the values i = x,y,z). We have 
decomposed the tensor structure of Eq. (10) into 
a vector and a spin part. This approach is now 
generalized to other corrections SH to the leading 
Zeeman Hamiltonian Hmo, 



(11) 



where 7 counts the correction terms. The opera- 
tors G 7 are linearly coupled to the magnetic field 
B. We split each element G 7 into a tensor struc- 
ture of spatial coordinates coupled to the magnetic 
field B as well as spin matrices a a . From the spin- 
less terms of the form v ■ B, we obtain the vector 
coordinates v l . 

The second-order spatial tensors A %3 in terms of 
the form A %1 u % a B 3 can be tensorially decomposed 
into a diagonal (scalar) part d, a symmetric tensor 
part t, and an antisymmetric tensor part r, 



A 13 = - S ij + t ij + r ij 



A 



A7v 



(12a) 
(12b) 

Is* A kk , (12c) 
S , (12d) 



where the summation convention is used for the 
Cartesian coordinates as is done throughout the 
paper. 

The orbital angular momentum part in leading 
order is identified as the vector term v a , and the 
spin part as related to the scalar operator do. For 
the finite mass correction in Eq. (5) with G\ in 
tensor form, we have the identification 



G '= - 5 E "L 

«l,a6 = ( f a X PbY 



b > 



(13a) 
(13b) 
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which is included as the first term (7 = 1) in SHm- 
This grouping is extended to higher-order terms 
and to make contact with the literature. According 
to Appendix A of Ref. [9] and Eqs. (2)— (4) of 
Ref. [4], we can split the c/l and gs factors into 
leading-order terms, denoted by the same symbols, 
and corrections SgL and Sg$, which, when added 
to gL and gs, yield the complete results g' L and g' s , 
which include the correction terms. So, for triplet 
P states, 

$9l = g' L - 9l , Sg s = g's - 9s ■ (14) 

This notation has been introduced in the theo- 
retical analysis of the experimental data for he- 
lium 2 3 P-states [9] based on angular momentum 
methods [4, 12, 13]. Compared to Eqs. (22) 
and (23) of Ref. [4], the prefactors in the expres- 



sions SgL = g'i — \J (2L + 1)L(L + l)/6 gh and 
Sgs = g' s — sj (25 + 1)5(5 + l)/6 gs evaluate to 
unity for triplet P states; for singlet P states, 
the first equality in Eq. (14) remains valid while 
9s — 0- The symmetric tensor parts t IJ are re- 
lated to the g x factor [9], and the mean values of 
the antisymmetric part r 13 result in a zero cor- 
rection. Later, these quantities were determined 
in subsequent theoretical calculations of other au- 
thors [4, 14]. We follow these conventions in order 
to be able to compare our final formulas with their 
results. We here use a Cartesian decomposition 
of the higher-order Zeeman Hamiltonian, as an al- 
ternative to angular algebra methods with 3 J, 6 J, 
and 9 J symbols [9], and identify the tensor contri- 
butions to <7l, gs, gx as described in the following. 



I 

C. Relativistic corrections 



Rclativistic corrections have been derived from the Breit Hamiltonian [15-17]. We follow formulas from 
Eq. (32) of Ref. [5] with six relativistic corrections to the Zeeman effect, 

5H Te i = hbol } j < - -y |_(r a xp a ) + a a \ ■ B + Z — -3 — (p a ■ a a ) (p a ■ B) > 

a a ^ 

, „ „2vl .95 - 1 (r ab x a a )(f a x B) g s (f ab x a h )(r a x B) p l a ( & , r l ab rj b \ g A 

+/iBQ k t — ^ — + ^w b + —) irbXB) \' 

SH iel = hbG 2 -B. (15) 
It is straightforward to identify the Cartesian tensor form of the relativistic correction G 2 , 

n i 2 f 1 f 1 , , is, z (gs ~ 1) f2 i Aj \ g s - 2 / d 2 , a < y A \ 
G 2 = f i B a 2^ \ 2 ( V ' 2 ' a ' a ° a ' 4 V3 3,a Ga ~ 3 ' a a ) 2~ \ — 3 — 4 ' a a j) 



1 / , 

2 



(16) 



Indeed, the tensor components from the first four 
terms in Eq. (16) read as follows, 



d 2 ,a 


-fa, 


(17a) 


Via 


= Pa ^a X PaT , 


(17b) 


d3,a 


1 

r a 


(17c) 


t ij 


= — ( r V' - -S ij r 2 1 
r 3 ^r a r a 3 r a j , 


(17d) 


f ij 
I 4,a 




(17e) 



Furthermore, we have the following terms from the 



fifth corrections in Eq. (16), 

d$,ab = Ta [ ab , (17f) 

r ab 

.ij _ 1 r a r Qb + r a r ab _ ^a [ ^ab Hj \ 

5 > ab ~rl b { 2 3 *)> 

(17g) 

riU = ^T (riH-rirl) ■ (17h) 
Finally, the sixth term in Eq. (16) yields a remain- 
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ing vector structure, 



III. EVALUATION OF THE g FACTOR 



61, ab 



I', 



62, ab 



(r a x p b y 

Tab 

{r a x r b ) 1 (r ab ■ p b ) 



(17i) 
(17j) 



We proceed to a final numerical evaluation of these 
corrections later. 



D. Self-energy correction 



In first-order perturbation theory, one requires 
only the diagonal matrix elements in the total an- 
gular momentum J. Insofar as first-order theory 
is concerned, one may replace Hmo + SHm by its 
restriction to the (2 J + l)-dimensional subspace 
spanned by the orthonormal vectors for to,/ = 
— J, — J + 1, . . . , J. Then, the magnetic Hamilto- 
nian linear in the B field can be rewritten as fol- 
lows, in terms of the gj-factor and J, 



We follow Rcf. [5] and base the calculation of 
the low-energy part of the self energy proceeds on 
a nonrclativistic Hamiltonian in the presence of an 
electromagnetic field in the length gauge, 



H = Hn + H 



A/0 



H~, 



efi-E-er 2 -E, (18) 



where Hq is the unperturbed Hamiltonian of the 
atom, H 1 is the Hamiltonian of the photon field, 
Hmo is the leading-order magnetic interaction 
given in Eq. (10), and the two dipole interaction 
operators describe the interaction of the bound 
electrons with the quantized electromagnetic field. 
The self energy has the form 

SE = - — / duuj 3 U\(r 1 +f 2 ) 
3tt Jo 

rr T rr 1 p — (n+f 2 )|0). (19) 

It is understood that 5 E is to be expanded in first 
order in the magnetic field B. Then, replacing the 
coordinates by electron momenta and using com- 
mutation relations, it is easy to rederive Eq. (38) 
of Rcf. [5] , additionally assuming that the state (f> 
has definite itll and ms quantum numbers (pro- 
jections of the orbital and spin angular momenta 
onto the quantization axis). We might just as well 
assume that the reference state has a defined value 
of the magnetic quantum number mj of the total 
angular momentum, 



SE = - 2 fi B — \ i drs ( 
7T I 



mj\(pi +P2Y 



y.\n\2{H -E)\{p 1 +p 2 ) s \m,,) \B 



(20) 

Adjusting the self energy correction to our conven- 
tion, we obtain the tensor structure 



SE = - 2/i B — (mj\G 3 ■ B\mj) 



(21) 



where G\ = v\ is given by the expression in curly 
brackets in Eq. (20). Thus, the self-energy correc- 
tion has a simple vector structure, and contributes 
to the orbital momentum L part. 



H 



MO ' 



5H 



= fi B ■ B = hb gj J- B . (22) 

7=0 



In order to calculate the gj factor, it is helpful to 
write it as an average over all magnetic projections. 
Using the shorthand notation \mj) for the state 
with quantum numbers \L, S, J, TUj), an important 
relation is 



(mj\G ■ B\mj) = Imj 



GJJB 



■J 2 



mj 



9J 



raj gj B , 

(mj\G ■ J \mj) 
J(J + 1) 



(23) 



This relation holds for any mj, and G may stand 
for any of the G 7 or for the sum G = X) 7 =o^7- 
Summing over mj and dividing by the number of 
states 2 J + 1, the gj factor can be determined as 
follows, 



1 



9J 



(24) 



J(J+1)(2J + 1) 
J 

(mj\Gi- J\mj) . 

i—0 mj=—J 

This expression involves a sum over the angular 
momentum projections and is manifestly indepen- 
dent of TO,/. 



A. Reduction of spin degrees and J 

The G 1 with 7 = 0, 1, 2, 3 have been defined in 
Sees. II B, IIC and II D. According to Eq. (12), the 
terms can be decomposed into diagonal d terms, 
and symmetric as well as antisymmetric t and r 
terms. Using Eq. (24), one can express the g- factor 
contribution of a term of the form G ■ J in terms of 
the multiplicative factor J(J+ 1)(2J+ 1) a correc- 
tion to the gj factor, and radial matrix elements. 



G 



For a contribution to G 7 • J of the form da ■ J, one 
can deduce for P states the formula 
J 

(m,j\ da ■ J \mj} = (25a) 

mj=—J 

A JS J(J+l)(2J+l) (^ k \d\^ k ), 

a 

A 21 = A 11 = -, A w = 0. (25b) 

The Cartesian basis \ip k ) of P states is normalized 
to (V> fe |^) = | 5 ki [see Eq. (31)]. The basis of the 
\tjj k ) states contains states without an explicit spin 
wave function, where the coordinate part is sym- 
metrized or antisymmetrized, according to Eq. (31) 
below. For a vector v coupled to J in Eq. (24), we 
use the following reduction scheme 

(mj\v ■ J\mj) = (26a) 

mj=—J 

Bjs J( J + 1) (2 J + 1) i e ijk (V>V 1^), 
B 21 =B 11 = ~, B w = ~. (26b) 

For a symmetric, traceless (quadrupole) tensor, we 
can project onto the Cartesian basis for P states 
as follows: 
,/ 

^ (mj|^V J J '|mj) = (27a) 

mj=—J 

Cj S J(J + i)(2j + i) J2 (i> j \t jk \^ k ), 

a 

C 2i = -^, G ^ = \ Gi O = 0. (27b) 

Finally, for a antisymmetric tensor r lJ coupled to 
a l a J 3 , the total contribution to the g j factor van- 
ishes for all states under investigation here. 

For excited helium P states, the leading-order 
expression (2) evaluates to 

gj(n 3 P J=0th2 ) = I , gj^Pi) - 1 , (28) 

where n is the principal quantum number of the ex- 
cited state, and the result for n 3 P is not of phys- 
ical interest because the magnetic projection for 
the state with J = always is fx = 0. The cor- 
rection Sgj to the Lande g factor can be expressed 
in terms of Sg^, 5gs, and g x [sec Eq. (14)], and 
prefactors Ajg, -B/s, and Gj$, 

Sgj= Aj S Sgs-2Bj S Sg L + ^Cj S Sg x , (29) 

where the Aj$, Bjs and Cjs coefficients are given 
in Eqs. (25), (26) and (27), respectively. It is prob- 
ably useful to note that our scheme is easily gen- 
eralized to other low orbital momentum states, for 



example, D-states, which have submanifolds with 
S = and S= 1, and J = 1,2,3. 



IV. NUMERICAL EVALUATION 

The nonrclativistic wave function of the P 
state ip and its energy Eq are determined for a 
Schrodinger equation with a nonrelativistic Hamil- 
tonian Eq. (8) 

H 4> = E iP, (30) 

based on the Rayleigh-Ritz variational principle. 
We use a basis set of explicitly exponentially cor- 
related functions (following Refs. [20, 21]) 

TV 
m— 1 

(31) 

which for the singlet (triplet) states is symmetric 
(antisymmetric) under an exchange of spatial coor- 
dinates, as required by the Pauli exclusion princi- 
ple. We rccmphasize that the superscript k denotes 
the Cartesian coordinate, that is, the wave func- 
tion with an orientation along the x axis would be 
denoted as %p x and involve the x coordinates rf and 
(in the exchange term) rf . The parameters (a, b, c) 
for the ith function are randomly generated from 
an optimized box (Ai,A 2 ) x (Bi,B 2 ) x (Ci,C 2 ) 
under the additional constraints a& + bk > s as 
well as 6fc + Cfc > e and Ck + ak > £, where 

e = 1J2 (Eq — E ) with Eq being the lowest sin- 
glet (triplet) P state energy for He + . 

In order to obtain a more accurate representa- 
tion of the wave function, we use two boxes that 
model the short-range and medium-range asymp- 
totics of the helium wave functions. In this basis, 
the matrix element of the nonrclativistic Hamilto- 
nian Hq can be represented as a linear combination 
of the integrals 

T(a, 6, C, 771,772,7712) (32) 
= J d'Vi d 3 r 2 rr 1 " 1 ^ 3-1 ^ 13-1 c -°ri-br 2 -cr 12 j 

with non- negative n\, n%, and n\2- Methods for 
their computation are well known [22]. The linear 
coefficients d m in Eq. (31) are obtained from a 
solution of a generalized eigenvalue problem. The 
numerical accuracy of the results is estimated from 
the apparent numerical convergence of the matrix 
elements as the size of the basis is increased. For 
the calculation of ip , we use an expansion with a 
moderate number up to 2 vV = 900 basis functions 
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TABLE I: Mean values of the tensor structures entering Eqs. (25), (26) and (27), for 2 1 P, 2 3 P and 3 3 P states. 
In view of Aio = do = 0, only the phenomenologically relevant results are indicated for singlet P states. 



Operator 


2 1 P 


2 3 P 


o 3 p 


ie !3fe (V'>o,alV' fc ) 


— 2.0 


— 2.0 


—2.0 


,-0.805 549 556(6) 


-1.096 171 714(2) 


-1.366 172(4) 


i Eak (Siblnn \8ibt r ,J) 


0.402 774 778(8) 


0.548 085 8571(4) 


0.683 086 2(2) 




-0.131044 018 6(5) 


0.256 875 920 7(3) 


0.069 756 861(4) 




-1.131383 6(3) 


-1.204 2321(4) 


-1.025 36(6) 


ieijki^lvi a \ip k ) 

(i> k \d 3 ,a\i> k ) 


A 10 = 


4.110 292 724 2 


4.116162 168 5 


-0.483 020 29131 


-0.216 764422 66 


-0.259 969 064 25 


A w = 


2.109 944 7016 


2.116 045 575 2 


Cio = 


0.072 236 399(2) 


0.080111516 9(4) 


Cio =0 


0.077 581379(2) 


0.084 694 997 4(3) 


& k \d 5 ,aty k ) 

W\ti k a b \^) 
ieiife(i|w6i, a JV fe > 


A w = 


0.109 596 679 06 


0.115 928 98192 


Cio =0 


0.066 891418 855(5) 


0.075 528 036 275(1) 


-0.077 009 223 65(3) 


-0.024 049 685(6) 


0.030 684 7511(4) 


0.029 866 744 798(1) 


0.008 854 469 40(2) 


-0.001 716 909 955(6) 




0.195 754(2) 


0.264 705(2) 


0.088 415(4) 



TABLE II: Sgj contributions to the singlet 2 Pi state 
where Sgs = 5g x = 0. The fine-structure constant is 
a = 1/137.035 999 074(44), and the electron-a mass 
ratio is m/M = 1.370 933 555 78(55) (see Ref. [18]). 
Theoretical uncertainties come from our estimate of 
higher order effects in the order a 4 (first parentheses), 
and finite-mass relativistic correction of order a 2 X (sec- 
ond). 



2 1 P 1 



Sg L x 10 6 



Finite mass [Eq. (35)] 
Relativistic [Eq. (15)] 
Self-energy [Eq. (21)] 
Total 

Theory: Ref. [3] 
Theory: Ref. [4] 



-8.968 94 
-7.85319 
0.024 22 
-16.798(9)(7) 

-15.771 
-16.810165. . 



"There is no uncertainty estimate given in Ref. [4]. 

(we use a prefactor 2 in order to clarify the distri- 
bution of the basis functions onto two variational 
boxes, as described in Ref. [21]). The numerical 
accuracy of the following nonrelativistic reference 
values, 

£0(2^)= - 2.123 843 086 4981013(3), (33a) 
E {2 3 P)= - 2.133164190 779 2831(2), (33b) 
£ (3 3 P)= - 2.058 081084 274 275(1), (33c) 

is fully sufficient for our calculations. All entries 
in Eq. (33) are consistent with the values given in 
Table III of Ref. [23]. For 2 3 P, our result also is in 
agreement with the entry in Eq. (20) of Ref. [24]. 

Using the wave functions ip k , we are able to ob- 
tain all necessary mean values defined by the tensor 



components d, v l and t lJ . In order to perform the 
calculation, the set of integrals of the form given 
in Eq. (33a) needs to be extended by additional 
classes including one or two of the indices ni, 122, 
and rii2 being equal to —1. The analytic formu- 
las for the master integrals and related recurrence 
schemes are well known, as given in Ref. [25]. Nu- 
merical values for the states with definite orbital 
momentum and spin-angular symmetries 2 : P, 2 3 P 
and 3 3 P are presented in Table I. Because the 
coefficients Aio an d Cio vanish, we do not pro- 
vide values for the d and t % i elements for singlet 
P states. Values without an error estimation are 
cut to eleven digits; all of these are believed to be 
numerically significant. 

In order to determine the finite mass effect of 
order 0(A 2 ), we first calculate the mass polariza- 
tion correction to the wave function, scaling the A 
parameter out of the perturbation, 



1 



(Eq — Ho) 



Pi -P2 



i> k ) ■ (34) 



The expression ISip^p) is relevant for the entries in 
the third column of Table I. The operator pi ■ P2 
in Eq. (9) changes neither the orbital angular mo- 
mentum nor the spin symmetry when acting on ip k 
in Eq (34). Thus, it can be expressed using a basis 
consisting only of the ip k defined in Eq. (31). Vari- 
ational parameters for Sip mp are generated in anal- 
ogy to those for the wave function tp, but the size 
of the basis is chosen to be larger (2N mp = 3 AT). 
With these results in hand, it is straightforward 
to calculate the mass polarization correction for a 
given operator. The only effect is the second-order 
correction to the v\ ab in Eq. (13). Together with 
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TABLE III: Sgj contributions to triplet P states. Again, the theoretical uncertainty of the final theoretical 
prediction comes from the estimate of higher order effects in the order a 4 , and finite- mass relativistic correction 
of order a 2 A (first and second parentheses, respectively). 



2 3 Pj 


5g L x 10 6 


Sg s x 10 6 


8g x x 10 6 


Finite mass [Eq. (35)] 
Relativistic [Eq. (16)] 
Self-energy [Eq. (21)] 
Total 


17.620 32 
-6.912 92 
0.032 74 
10.740(2)(2) 


0.0 
-80.429 3 
0.0 

-80.43(2)(4) 


0.0 

-5.385 487 
0.0 

-5.385 5(13)(2) 


Theory: Ref. [4] 
Theory: Ref. [2] 
Theory: Ref. [3] 
Theory: Ref. [5] 
Experiment: Ref. [1] 
Experiment: Ref. [9] 


10.719 291 ..." 
10.6(4) 

8.838 
10.752 033 . . . 

4.9(1.9) 

3.8(9.0) 


-80.436 904. . . a 

-80.46(1) 

-80.401 

-76.0(2.4) 


-5.391 808... a 

-3.5(1.5) 

-5.344 

4.0(25.0) 


3 3 Pj 


5g L x 10 6 


8g s x 10 6 


Sg x x 10 6 


finite mass [Eq. (35)] 
relativistic [Eq. (16)] 
self energy [Eq. (21)] 
Total 


4.788 03 
-3.029 57 
0.010 94 
1.769(9)(2) 


0.0 
-75.083 9 
0.0 

-75.08(2)(4) 


0.0 
-2.648 665 
0.0 

-2.648 7(7) (3) 


theory: Ref. [4] 
theory: Ref. [19] 


1.772 223...° 
-0.17(2.8) 


-75.096 557...° 
-75.13(3.27) 


-2.650 192..." 
-2.75(10.02) 



"There is no uncertainty estimate given in Ref. [4]. 



part of the third order correction to Vq a , we obtain 



<<mi^) = --£ [m4, ab \^ k ) (35) 

a 

where the formula is expressed in compact form by 
writing the coefficient of the first term as m/M, 
not n/M . The other terms in the third order per- 
turbation of Vq a result in the shifted g^ coefficient 
in the leading term in Eq. (5) 

= 0.999 862 916 942 649(5)(55), (36) 

where the first uncertainty estimate refers to the 
numerical uncertainty of the particular contribu- 
tion (finite-mass correction to gi/), and the second 
uncertainty comes from the CODATA electron-a 
mass ratio m/M = 1.370 933 555 78(55). 

The most numerically intensive part of the g fac- 
tor calculations in Eq. (24) is the vector component 
of the self-energy correction (21), which can be re- 
duced [5] to the expression 

ky^Vr |^) = (S ir 5 ks - S is 5 kr ) <^|(pi +AY 
x\n\2(H Q -E)\(p 1+ p 2 ) s |^), 

(37) 



which differs from the ordinary Bethe logarithm 
in the absence of a linear term Hq — E multiply- 
ing the logarithm. Matrix elements involving the 
logarithm of the Hamiltonian necessitate the use of 
the methods usually employed for Bethe logarithm 
calculations for excited states in helium, where due 
to conceivable numerical challenges, for a long time 
asymptotic formulas [26-29] in l/n and 1/Z were 
the preferred method of calculation. Direct and 
accurate calculations of logarithmic sums over the 
helium spectrum have become possible only quite 
recently [30-32]. Here, we closely follow to the in- 
tegral representation of the Bethe logarithm [33], 
which for the expression in Eq. (37) has a particu- 
lar compact form, 



(38a) 



t 3 ' 

/(f) = (6 ir 5 ks - 6 iS 8 kr )(ij l | (pj + p 2 ) r (38b) 

(pi+p 2 ) s |^ fc ), 



1 



t 



Ho-E +uj 
1 



(38c) 



We perform an integration over 100 equally spaced 
and optimized t points, following ideas outlined in 
Ref. [33]. The well-defined limit of the integrand, 
lim^o f(t)/t 3 = 0, facilitates the numerical eval- 
uation. 
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V. RESULTS 

The numerical data for the individual operators 
allows us to obtain theoretical analysis of the L, S 
and the x part of the Lande g factor, as defined in 
Refs. [4, 9] and discussed above. We express our 
results in terms of Sgj, which is obtained as the 
difference of the total prediction and the leading- 
term [see Eq. (14)]. For the numerical evaluation, 
we use Eq. (24). We express the correction 5gj 
to the Lande g factor in terms of SgL, Sgs, and 
Sg x , and prefactors Ajs, Bjs, and Cjs, as given 
in Eqs. (25), (26) and (27), and (29). 

We keep the conventions of Refs. [4, 9] and com- 
pare our results to the experimental and theoreti- 
cal literature. In Tables II and III, we provide data 
split into a finite-mass part related to Eq. (13) in- 
cluding the second order mass polarization correc- 
tion Eq. (35), the relativistic correction Eq. (16), 
and the self-energy correction given in Eq. (21). 
The first conceivable source of uncertainty for these 
contributions is purely numerical, due to the finite 
numerical accuracy of the components in Table I. 
However, in most cases, the numerical uncertainty 
is negligible as compared to the omitted higher- 
order effects. 

We are able to report that our results confirm 
the numerical data reported previously in Ref. [4] 
at the level of relativistic operators without finite- 
mass corrections. Here, we attempt to go beyond 
the leading relativistic effects. The self-energy cor- 
rection to the P state g factor consists of two parts, 
one of which involves a Bcthc-logarithm type term 
(logarithmic sum over virtual excited states) and 
is ultraviolet finite, in contrast to the Bethe log- 
arithm contribution to the S state Lamb shift, 
which is known to be ultraviolet divergent [34, 35]. 
The second contribution due to the self-energy is 
a high-energy contribution, which is manifest in 
the anomalous magnetic effects of the free elec- 
tron, which are included into the relativistic Zee- 
man Hamiltonian given in Eq. (15). These con- 
tributions are infrared finite. We obtain very 
good agreement with a numerical result reported 
in Ref. [5] for the self energy correction Sg^ for 
2 3 P. Combining relativistic and radiative effects, 
we should mention the presence of an additional 
third term with a prefactor gs — 2 in Eq. (15). In 
Ref. [4], this term had not been taken into account, 
apparently, but its numerical magnitude does not 
shift the final result significantly. Our self-energy 
correction, which we add to the relativistic result 
of Ref. [4] , includes the relativistic anomalous mag- 
netic moment effects of the electron and the spin- 
dependence of the Bethe logarithm term. We also 
include 0(A 2 ) corrections to the leading order re- 
sulting in the finite mass correction in Eq. (35). 
Such correction gives —0.8 x 10~ 8 , —1.0 x 10~ 8 



and -1.3 x 10~ 8 to the Sg L of 2 1 P, 2 3 P and 3 3 P 
respectively, having the order of the self energy 
correction. 

There are two significant sources of theoretical 
uncertainty for the final results, which are given 
by finite-mass corrections to the relativistic effects 
and by higher-order QED contributions. In the 
results reported in Ref. [4], the mass scaling and 
mass polarization corrections to the relativistic ef- 
fects have been included; these results are of rel- 
ative order 0(a 2 A) with respect to the leading g 
factor term. However, there are additional finite 
mass relativistic effects of the same order, which 
can be deduced from Eq. (40) in Ref. [11], and 
which should be included in a systematic treat- 
ment. Here, we do not perform a complete cal- 
culation of the terms of order 0{a 2 A), and so we 
do not include the relativistic reduced-mass correc- 
tion at all in our final results. We use some partial 
results we have obtained in the order 0(a 2 A) in 
order to estimate the size of the relativistic-recoil 
correction. 

These include the scaling and the mass polariza- 
tion corrections for the L part in 2 1 P\ (4.4 x 10 -9 ), 
as well as for the x part (1.6 x 10~ 9 ) and the S 
part (2.5 x 10~ 8 ). For 2 3 P 7 , we have 1.1 x 10~ 10 
for the x part, 9.2 x 10~ 10 for the L part, and 
2.5 x 10 -8 for the S part. Finally, for the x part in 
3 3 Pj, we have a result of 1.7 x 10 -10 . These results 
guide our estimates of the theoretical uncertainty 
indicated in Tables II and III, where we multiply 
the partial results with a conservative weight factor 
of 1.5. For the higher-order QED contributions, we 
are not even able to present approximate formulas 
based on, for example, hydrogenic contributions, 
because the theory of the Zceman effect has not 
been developed until now to this order. Therefore, 
we use the combined value from the relativistic cor- 
rection multiplied by a factor a 2 , and the leading 
QED correction multiplied by a factor a, to obtain 
a conservative estimate of the uncertainty due to 
the uncalculated higher-order effects. We employ 
an additional conservative enlargement factor of 5 
in order to estimate the size of the effects of or- 
der a . 

As evident from Tables II and III, agreement of 
theory and experiment is rather satisfactory for a 
number of contributions, with the exception of a 
2 a discrepancy for Sgs (comparing our result to 
that of Ref. [9]) for the 2 3 P state and a 3cr dis- 
crepancy for SgL (comparing our result to that of 
Ref. [1]). It would be very interesting to remeasure 
the effect and clarify the status of the experimental 
results. 
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VI. CONCLUSIONS 

In the calculation of the gj factor of excited 
states of helium, all aspects of atomic physics play 
a role: electron correlation, relativity, and QED 
radiative corrections. Furthermore, these effects 
are all intertwined, that is, there are QED radia- 
tive corrections multiplying the relativistic effects, 
contributing in higher order. We here carry out a 
theoretical analysis of the P state g factor of singlet 
and triplet helium states, with a special empha- 
sis on relativistic, radiative, and finite-mass cor- 
rections. The effects arc calculated through order 
a 2 (for the relativistic effects), and we also include 
radiative effects of order a 3 . Furthermore, finite- 
mass corrections of order A and A 2 are included, 
while available partial results for the effects of rel- 
ative order a A are used in order to estimate the 
theoretical uncertainty in this order. Results are 
summarized in Tables II and III; a 2 a discrepancy 
for Sgs in comparison to Rcf. [9] and a by now- 
famous 3<T discrepancy for 5gL in comparison to 
Ref. [1] highlight the need for additional experi- 
mental evidence before definitive conclusions can 
be drawn. 

With regard to the QED self-energy correction 
to the g factor, one encounters a peculiar situation 
for P states: Namely, both the high-energy part 
as well as the low-energy part arc separately fi- 
nite, while the low-energy part is given by a Bethe- 
logarithm term, and the high-energy part is given 
by the anomalous magnetic moment. Both terms 
are finite, and it is easy to overlook one of the con- 
tributions. For atomic hydrogen, this has been ver- 
ified both numerically [7] and analytically [8]. For 
helium singlet versus triplet states, we find that 
the low-energy part of the self-energy correction 
to the bound-state g factor is spin-dependent. 

Our calculations are performed in an angular- 
momentum coupling scheme which allows us to 
separate the internal degrees of freedom of the 
atom from the interaction with the external mag- 
netic field. The angular momentum algebra can 
become rather involved for helium P states. We 
fully confirm the relativistic treatment of Ref. [4] 
using our mixed approach. In a more general con- 



text, one may recall that the contribution of the 
anomalous magnetic moment to the P state life- 
time in few-electron systems has recently given rise 
to interesting experimental-theoretical discrepan- 
cies [36] which remain to be fully resolved. 
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Appendix A: Nonrelativistic Treatment 

The nonrelativistic Hamiltonian of an n-electron 
atom is given as (in atomic units) 



H 



Pn 
2M 



E 



2m r a 



E— ' ( A1 ) 



a>b 



where we keep the electron mass m and the nucleus 
mass M in symbolic form. In the center-of-mass 
system, we have pn = ~^2 a Pai and therefore 



(A2) 

where the latter term corresponds to the mass po- 
larization. The reduced mass fi is given as 



1 , Pa-Pb 



1 _ 1 1 

fx m M 



(A3) 



If we define the ratio A = —/j,/M as in Eq. (4), 
then an important identity is 



M m 



(A4) 
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